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Using ultrafast optical techniques, we detect two types of bosons strongly coupled to electrons in 
the family of Bi2Sr2CaCu208+5 from the underdoped to overdoped regimes. The different doping 
dependences of the electron-boson coupling strengths enables us to identify them as phonons and 
spin fluctuations: electron-phonon coupling (Ae-ph) peaks at optimal doping, and electron-spin 
fluctuation coupling (Ae_s/) decreases monotonically with doping. This observation is consistent 
with two facts: (1) superconductivity is in close proximity with antiferromagnetism at low dopings, 
and (2) a pronounced lattice renormalization effect at larger dopings. 

PACS numbers: 



Despite many advances in understanding copper-oxide 
high-transition-temperature [Tc) superconductors, there 
still exists no universally accepted mechanism. De- 
termining the nature of interaction responsible for the 
Cooper-pair formation remains one of the grand chal- 
lenges in modern condensed matter physics. The most 
probable candidates are lattice vibrations (phonons) 
[U [2] , spin fluctuation (SF) modes [3] , and pairing with- 
out invoking glue [J. For conventional superconductors, 
structure in the electron tunneling dl / dV characteristics 
established unambiguously that the attractive pairing in- 
teraction was mediated by phonons [5 . For high- Tc su- 
perconductors, structure in dl /dV has also been found 
in many tunneling measurements [5]. More recent scan- 
ning tunneling microscopy (STM) experiments revealed 
an oxygen lattice vibration mode whose energy is anti- 
correlated with the local gap value on hole-doped Bi-2212 
[7] while a bosonic mode of electronic origin was found 
in the electron-doped Pro.88LaCeo.i2Cu04 [8]. Together 
with salient features observed in angle-resolved photoe- 
mission spectroscopy ( ARPES) [2^ ^ , these new results 
raise the fundamental question of whether the bosonic 
modes are a pairing glue [lOj or a signature of an inelas- 
tic tunneling channel }TT| . 

The role of the electron-boson interaction in high- 
Tc superconductors has been studied by different tech- 
niques. For example, inelastic neutron scattering tracks 
the changes in boson energies or dispersions upon en- 
tering the superconducting state. ARPES \TT measure 
the effects of electron-boson interaction on electronic self- 
energies, and planar junction experiments determine the 
energy of the bosonic mode [T3] . STM measure the local 
density of states through the local differential tunneling 
conductance, where the characteristic boson mode energy 
is estimated from the dip position f7^. However, it does 



not give the electron-boson coupling strength directly be- 
cause both the coupling strength and mode energies are 
encoded in the electron self-energy itself. 

As for the electron-boson coupling strength in 
cuprates, time-integrated optical measurements can be 
useful [m [15] , but it is difficult to elucidate whether one 
or more bosons are involved. Time-resolved pump-probe 
spectroscopy is a powerful technique used to probe the re- 
laxation dynamics of photoexcited quasiparticles in cor- 
related electron systems such as cuprate, pnictide and ac- 
tinide superconductors, spin density wave materials, and 
Kondo systems [TM2^ . Its unique contribution lies in 
its ability to extract the value of the electron-boson cou- 
pling strength (A) directly, via the electron-boson relax- 
ation time, without the need to perform complicated in- 
version algorithms. This procedure has been experimen- 
tally verified on the conventional superconductors [23) . 
In this Letter, we report measurements of time-resolved 
quasiparticle relaxation of high-quality single crystals of 
underdoped (UD) to overdoped (OD) Bi2Sr2CaCu208+5 
(Bi-2212, hole concentration p=0. 10-0.22). Our data in- 
dicate the coupling of electrons to two bosonic modes: 
the electron-phonon coupling constant (Ae-p/t) peaks at 
optimal doping, while the electron-SF coupling constant 
(Ae_s/) decreases monotonically with doping. 

The family of the bi-layer cuprate Bi-2212 has been 
the most intensively studied class of high-T^ supercon- 
ductors in recent years, due to their (a) extreme cleav- 
ability, (b) containing only Cu02 planes and not chains, 
and (c) the possibility of growing samples with a larger 
range of T^s. Single crystals of Bi-2212 were obtained 
from three groups (two Tokyo and one Tsukuba) grown 
by the floating zone method with doping control, yield- 
ing values of T^ (determined by magnetization data) that 
depend on the hole doping level (p) spanning from the 
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FIG. 1: (a) AR/R versus pump-probe time delay, of an UD 
Bi-2212 sample {Tc=65 K). (o): Experimental data. Blue 
line: one-exponential fit. Red line: two-exponential fit. (b) 
Fast relaxation rate l/rfast versus Tc- Note the difference in 
the strength and sign of the correlation in the OPT-OD and 
UD samples. Red dotted line is a guide to the eye. (c) l/rfast 
versus doping. In this representation l/r/ast decreases mono- 
tonically with increasing doping, but with a faster decrease in 
the UD region, (d) 1/tsIow versus doping. 



UD (p=0.10, Tc=65 K) to the OD (p=0.22, Tc=65 K) 
regime. In the optimally-doped (OPT) sample, Ca has 
been doped with Y to obtain the highest Tc of 95 K. 
Underdoping was achieved using excess Bi atoms substi- 
tuted for the Sr sites as well as reducing oxygen content, 
while the more OD samples have been doped with Pb, 
to obtain lower values of T^. The values of p were ob- 
tained from the T^ values using the parabolic law [23] 
T^lT^ax ^ -^ _ 82.6(p _ 0.16)2, where r™°^=95 K. 

To avoid any competing relaxation processes from 
emergent low temperature states (e.g., superconducting, 
pseudogap, antiferromagnetic, or stripe order), we take 
all data at room temperature. Figure fit a) shows the 
time dependence of the photoinduced change in reflection 
(AR/R) of an UD Bi-2212 sample (see supplementary in- 
formation |25) for more discussion). The time evolution 
of the photoinduced reflectivity change AR/R first shows 
a rapid rise (of the order of the pump pulse duration) fol- 
lowed by a subsequent decay. As shown in Fig. llTa) , the 
data can be fitted better by two exponentials (red line) 
than a single exponential (blue line). It indicates the 
quasiparticle relaxation has two components — the fast 
component Tfast (~100 fs) and the slow component Tgio^ 
(~700 fs). We attribute the fast relaxation process to the 
electrons first transferring energy to a bosonic mode (e.g. 
phonons) which are more strongly coupled at a charac- 
teristic time Tfast- These bosons then continue to cool 
by energy dissipation via anharnionic decay at a charac- 
teristic time Tsiow 

Since the transfer of electron energy first occurs 
through selected modes that are most strongly cou- 
pled to electrons, we use Tfast to be indicative of the 



FIG. 2: Grey lines: Individual data sets. Open circles = 
averaged data, multiplied by a proportionality factor. The 
proportionality factor [(5.0 ±0.1) x 10^ K] is a fitting param- 
eter constrained to be doping independent, since the same 
amount of incident laser power is absorbed by all the sam- 
ples. Solid lines are time evolutions of the temperatures of 
the various subsystems: electrons (black), hot phonons (Tj 



red), cold phonons (T^^ , green), hot spin fluctuations {T!. 
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blue), and cold spin fluctuations {T^f'^, magenta). All data, 
taken at 295 K, are shown here to 3 ps, but were taken and 
fitted up to 7 ps. (d) Doping dependence of Ae_ph (blue) and 
\e-sf (red). Fittings to the individual data sets yield a dis- 
tribution of fitting parameters, and the standard deviations 
of these distributions give the error bars. 



electron-phonon coupling strength. These strongly cou- 
pled phonon modes should be the most relevant in dis- 
cussing the possible phonon-mediated superconductivity. 
Figure fl|b) shows the Tc dependence of the fast relax- 
ation rate 1/Tfast — there is a weak and positive corre- 
lation between 1/Tfast and Tc for OPT to OD samples, 
while for the UD samples, the correlation is strong and 
negative, in contrast to what is expected from the BCS 
theory, where Xe-ph (proportional to I /Tfast) correlates 
positively with Tc- Figure [l|c) shows the doping depen- 
dence of the fast relaxation rate 1/Tfast — here 1/Tfast 
decreases monotonically with increasing doping, with a 
faster decrease in the UD than the OD region. These two 
panels suggest that, in the UD region, another impor- 
tant relaxation mechanism, in addition to the electron- 
phonon interaction, contributes to the quasiparticle re- 
laxation. For the slow relaxation rate I/tsIow, as shown 
in Fig. fird), it varies only slightly with doping, but a 
little more pronounced at low dopings. 

We further elucidate the contribution of possible mul- 
tiple bosonic modes to our data and and extract their 
respective coupling strength to electron based on an ef- 
fective temperature model. The existence of multiple re- 
laxation rates in a single sample suggests a conventional 
two-temperature mode is not sufficient. We note that 
a three-temperature model (3TM) has been used in the 



time-resolved ARPES analysis of an OPT Bi-2212 sam- 
ple [26]. In the 3TM, photoexcited electrons first trans- 
fer their energy to the ("hot") phonons that are more 
strongly coupled with them. These hot phonons then 
lose their energy to the cold phonons through anhar- 
monic coupling. This model has been successfully used 
to fit the data for one single sample OPT Bi-2212 in 
Ref. 1261 However, to interpret the relaxation phenom- 
ena across the whole range of dopings as presented here, 
we need to go beyond the 3TM. Especially, the presence 
of strong SFs in the UD regime leads us to consider a 
five-temperature model (5TM). In this new model, there 
exists a second cooling channel for the hot electrons, that 
is, via coupling to hot SFs, which subsequently also cool 
via scattering with cold SFs. The time evolution of the 
temperatures of the various sub-systems: electronic (Tg), 
hot/cold phonon {T^hl^ph)^ hot/cold SFs {T'^f /T^^f) 
satisfy the rate equations: 

dT, _ P ^ nf ~ n^f 3Xsf{n'^fr n^f - n'^f 



dr 


Ce 


^' T. 




hnkg 


Te 


\^) 


dr C^ot 1 




nX 


rphot 


rpcold 

Tph 


(2) 


dT^f 


( 


Ce 3A,/(17 


hot\3 ^ 


'J- 


n'lf 


rj^hot _ y 


cold 


dr 




2 
B 


Te 




(3) 






dTX" 
dr 


/-'hot 
r<c.old 

'^ph 


rphot 

^ph 


rpcold 

^ph 
'Tph 




(4) 






dr 


r<hot 

(-icold 


rphot 


__ rpcold 




(5) 



The system is excited by a Gaussian pulse P with FWHM 
45 fs and energy density of 0.166 J/cm'^. The spe- 
cific heat of electrons, hot/cold phonons, and hot/cold 



SFs are Ce = lT,,C';'^' = Nphfph^XT^^^v"''' 
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(jcou ^ ^ (^ _ f )n-f90^^ respectively The pa- 

rameter /p^ ifsf) denotes the fraction of total phonon 
(SF) modes that are more strongly coupled to the elec- 
trons, and Nph {Ngf) denotes the number of phonon (SF) 
modes in the irradiated volume. The distribution func- 



tions are n 
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tions, but are results of performing delta-function energy 

integrals [27] . 

The energy of the hot phonon mode Slp^*=40 meV cor- 
responds to the out-of-plane out-of-phase oxygen buck- 
ling B\g phonon. Though cuprate samples like Bi-2212 



are inhomogeneous both in energy gap and characteristic 
boson frequency, the spatial average of mode frequency 
is doping independent [?]• Therefore, we assume ri\°* is 
constant throughout the entire doping regime. Neutron 
scattering data [28| reveal an accoustic phonon mode at 
20 meV. Since recent STM data I7| did not observe any 
coupling between electrons and this particular accoustic 
mode, we choose it to be the energy of our cold phonon 
bath (f^^^"). The hot (r^Jf) and cold {^1°^'^) SF energies 
are taken from optical spectroscopy [291 and are doping 
dependent: both il^,'?* and ^f}'^ scale with the centroid 
position of the broad bosonic background, with the con- 
straint that r2^^*=41 meV at OPT. Our model also takes 
into account the fact that (a) the samples change from 
a good metal at higher doping to a bad metal at lower 
dopings (via factor C\) [551 HO], and (b) the penetra- 
tion depth at 800 nm pulse changes with doping [29]. 
Two other parameters do not change with doping: (a) 
the Sommerfeld coefhcient 7 [5T1 W2\ and (b) fraction of 
incident laser power absorbed by the sample [2S. The 
two most important parameters to be extracted from this 
model are the coupling constants \e-ph and Ae_s/. 

Figure l2la)-(c) shows the STM fits of three repre- 
sentative samples. The time evolution of the various 
sub-systems (electron, hot/cold phonons/SFs) depends 
on the relative strengths of the coupling constants [25] . 
The optimal hot phonon and SFs were obtained using a 
parameter space study j^Sj. The change in peak shapes 
is well described by our model, and we obtained very 
good fits for the entire doping range. Figure [2Fd) shows 
the doping dependence of the respective coupling con- 
stants. In the UD region, \e-ph initially increases with 
doping between p=0.098 and 0.11, then starts to de- 
crease between 0.11 and 0.13 doping, before increasing 
again to a peak at OPT (p==0.16), and thereafter de- 
creasing with increasing doping in the OD region. On 
the other hand, \e-sf decreases with increasing doping, 
with a stronger decrease in the UD than the OD regions. 
We also performed fits using other possible hot phonon 
and SF modes: (a) half-breathing in-plane copper-oxygen 
bond stretching phonon as the hot phonon mode, with 
fip)J*=70 meV, (b) magnetic resonance mode as the hot 
SF mode, with doping-dependent energy values obtained 
from neutron scattering data [33j , and/or the peaks in the 
bosonic spectrum in optical spectroscopy data [2S1 [52] ■ 
Only the magnitudes of Xe-ph and Ae-s/ differ; the over- 
all trends with doping do not change. 

The monotonic decrease of Xe-sf with doping is con- 
sistent with the weakening of SFs with increasing dop- 
ing, resulting in a weaker coupling between electrons and 
SFs. Note the strong decrease of Xe^sf occurs up to 
at p=0.13; thereafter the decrease becomes more grad- 
ual. This change in behavior at p=0.13 could be related 
to the formation of stripe ordering [M] at p=l/8. The 
non-zero value of Ae-s/, even in our most OD sample 
(p=0.22), suggests that SFs are present even in the OD 
regime, in agreement with studies by inelastic neutron 
scattering [55] . 



For Xe-ph, its initial increase between 0.098 and 0.11 
doping is qualitatively consistent with BCS theory, where 
an increased electron-phonon coupling gives rise to an 
increased Tc and vice versa. Its subsequent decrease be- 
tween 0.11 and 0.13 (^1/8) doping, increase away beyond 
p^l/8, may be due to the presence of stripe order at this 
doping level |36) . The concurrent decrease of X^^sf in 
this same "stripe" region implies a suppression in the 
superconductivity, consistent with a plateau in the su- 
perfluid density near this 1/8 doping f37j. The doping 
dependence of Xe-ph in the OPT-OD region — maxi- 
mum at OPT, and decrease with overdoping, is also in 
qualitative agreement with BCS theory, and strong lat- 
tice renormalization effects [21 138] . 

Independent confirmation of the electrons coupling to 
multiple bosonic modes come from frequency-resolved 
and time-resolved pump-probe measurements of OPT Bi- 
2212 [39] . The resulting spectral function necessitates the 
electrons to be coupled directly to hot phonons, cold lat- 
tice, as well as a bosonic mode of electronic origin, which 
the authors suggested could be SFs or current loops. 
Moreover, ARPES data on heavily UD La2_a;Sr2.Cu04 
showed fine structure in the electron self-energy, demon- 
strating the involvement of multiple boson modes in the 
coupling with electrons [40]. The contribution of our 
work, besides working on a different class of cuprate 
superconductors, is that we have obtained the doping 
dependence of the different electron-boson coupling con- 
stants from the UD to the OD regimes. 

It is important to realize that our 5TM requires both 
the electron-phonon and electron-SF interactions to be 
responsible for the fast (~100 fs) relaxation. Impor- 
tant confirmation of this comes from the temperature- 
dependent pump-probe data of our most OD sample 



(Tc^65 K) '^. Besides observing the opening up of 
a pseudogap at T*~100 K, the fast relaxation time rjast 
increases with decreasing temperature, before peaking at 
T* and decreasing to ^100 fs at 30 K. This change in 
behavior of Tfast at T* is intriguing — it suggests that 
the electron-SF coupling, in addition to electron-phonon 
coupling, is involved in the initial fast relaxation of the 
hot electrons. The peak at T* , and its subsequent de- 
crease below T* , is then due to an increased scattering 
rate between electrons and SFs as the sample enters the 
pseudogap phase. This scenario is further confirmed by 
the temperature-dependence of Tfast above T* — a fit 
to l/T" yields n=1.3, which disagrees with the behav- 
ior predicted for the electron-phonon relaxation time for 
good {n = 2) and poor {n = 3) metals [50] . 

In conclusion, we have performed ultrafast optical 
techniques on a wide range of doping levels in Bi-2212. 
The relaxation analysis suggests the existence of two 
types of bosonic modes strongly coupled to electrons. 
The different doping dependences of the electron-boson 
coupling strengths enables us to identify them as phonons 
and spin fluctuations: electron-phonon coupling (Xe-ph) 
peaks at optimal doping, and electron-spin fluctuation 
coupling (Xe-sf) decreases monotonically with doping. 
The observation should shed new insight into the mech- 
anism of high-Tc superconductivity in cuprates. 
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Supplementary Information for 
Doping dependence of the electron-phonon and electron-spin fluctuation interactions in the 

high-Tc superconductor Bi2Sr2CaCu208+5 



In this EPAPS auxiliary online document, we provide (A) details of experimental pump-probe 
setup, (B) general comments on the temperature models, (C) a comparison between a five- and a 
three-temperature model, (D) doping-dependent and independent parameters in the 5TM, (E) 
curve-fitting procedure, (F) time evolution of the electronic, phononic, and spin fluctuation 
subsystems based on the fits, (G) details about the error analysis and (H) doping dependences of 
the electron-phonon and electron-spin fluctuation coupling constants using three other scenarios 
of the phonon and spin fluctuation mode energies. 

A. Femtosecond pump-probe setup 

In our experiments an 80-MHz repetition rate Ti:sapphire laser produces 45 -femtosecond 
pulses at approximately 800 nm (~1.5 eV) as the source of both pump and probe optical pulses. 
The pump and probe pulses were cross-polarized, with a pump spot diameter of -60 [xm and 
probe spot diameter of -30 |im. The reflected probe beam was focused onto an avalanche 
photodiode detector. The pump beam was modulated at 1 MHz with an acoustic-optical 
modulator to minimize noise. The experiments were performed with an average pump power of 3 
mW, giving a pump fluence of -0.1 |iJ/cm^ and a photoexcited quasiparticle density of 0.02/unit 
cell, showing that the system is in the weak perturbation limit. The probe intensity was -10 times 
lower. Resolution is at least 1 part in 10*". The fitted values of x have a typical error of ±1%. All 
data were taken at room temperature. For all the dopings in this paper, the same amount of the 
incident laser power is absorbed by the sample, i.e. -90%. 

B. Temperature Model 

On the timescale of several tens of femtoseconds (fs), it is not inappropriate to describe 
the various sub-systems with their corresponding temperatures. In the theoretical paper by Allen 
[1], he suggested a way to calculate the electron-phonon coupling constant (^e-ph) from the 
electron-phonon relaxation time (Xeph) in conventional superconductors. Brorson [2] followed up 
on his work by fitting pump-probe data of many conventional superconductors to the 
two-temperature model, and obtained values of the Ig^ph that agree very well with the literature. 
This shows that the effective temperature description of electrons is reasonable even in time 
scales of -100 fs. In the cuprate superconductors, the on-site Hubbard interaction is [/= 8?« 1.04 
eV, where t is the hopping amplitude [Refer to Ref. [3] to get f = 130.5 meV]. This corresponds to 
an electron-electron thermalization time le-e of -4 fs (using T~h/U , h being the Planck's 
constant). Therefore, owing to this large U, the electron-electron interaction is much stronger than 
that in conventional superconductors. This Te_e is much smaller than electron-boson long-tail 
relaxation time. As such, it will take a much shorter time for the photoexcited electrons to 
equilibrate to attain a Fermi-Dirac distribution. We maintain therefore, that it is reasonable to talk 
of an effective electron temperature in the cuprates at the sub- 100 fs timescale. Perfetti et al. also 
used the concept of electronic temperature in the interpretation of his time -resolved ARPES data 
[4]. Finally, the extracted electron-phonon relaxation timescale (sub-100 fs) itself makes the 
effective temperature internally consistent. 

C. Comparison: Five-Temperature (5TM) vs. Three-Temperature Model (3TM) 

Following Perfetti et al. we started out fitting our reflectivity data to a three-temperature 
model (3TM), where photo-excited electrons relax back into the ground state via electron-phonon 
interaction. In this section we explain how the 3TM failed to fit our data. 



When using the 3TM to describe multi-level electron relaxation processes, we assume that the 
measured reflectivity of the sample is proportional to the electronic temperature. This 
proportionality constant should be the same for all samples, and thus we do not treat it as a fitting 
parameter. Instead, we first tune this proportionality constant by hand for a sample near optimal 
doping, to pick a value that would enable us to get good fits. As it turned out, we can pick one 
value for this proportionality constant (pO, so that the tail of the electronic temperature-delay can 
be fitted out to 7 ps, at the expense of a good fit around the peak electronic temperature (see 
Figure 1). For another value of the proportionality constant (pa), we can obtain a good fit 
around the peak electronic temperature, but will no longer be able to fit the tail well out to 7 ps 
(see Figure 2). 
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Figure 1. (Nearly Optimal Doping). The fit shown in the left figure 
underestimates the peak. At the tail the three temperatures become equal i.e. the 
electronic system is in equilibrium with the cooling phonon bath (blue dotted line: 
experimental data multiplied by proportionality constant, black line: fitted 
electronic temperature, red line: hot phonon temperature, green line: cold phonon 
temperature). The right figure shows the absolute of the residuals at each point of 
the curve. The residuals become smaller at later times. 
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Figure 2. (Nearly Optimal Doping). The fit is in very good agreement around the 
peak (left). Equilibrium between the electronic and bosonic systems is reached 
around 2 ps. All temperatures remain constant from then on and thus cooling at 
later times is not captured. In comparison to Figure 1 we find that residuals (right) 
are now smaller around the peak but increase at later times. 



When we fit the ultrafast data of the other samples to the 3TM, using the proportionaUty constant 
Pi, we find the fits are poor both around the peak and the tail (see for example, the fit to an 
overdoped sample in Figure 3). The same is true if we use the proportionaUty constant p2 (see for 
example, the fit to the same overdoped sample in Figure 4). 
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Figure 3. (Overdoped Sample). If we use the same ratio between reflectivity and 
temperature as in Figure 1 for overdoped samples we find that the peak is not 
estimated properly. From the residuals (right) and from the actual fit (left) it can 
be seen that the tail is not properly captured as well. 
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Figure 4. (Overdoped Sample). If we use the same ratio between reflectivity and 
temperature as in Figure 2 we get an improved tail section but the peak still is 
strongly underestimated. 



In fact, no choice of proportionality constant will produce an acceptable compromise across the 
series of samples between fitting the peak electronic temperature, and also fitting the tail out to 7 
ps. This suggests that we need cooling channels beyond the one within the 3TM to fit the peak 
and the tail well across all dopings. We found that the 5TM, which contains two cooling channels, 
is the simplest extension of the 3TM that admits good fits to the peak and the tail across all 
dopings using a common proportionality constant. 



D. Doping-dependent and doping-independent parameters in the 5TM equations 

(1) Change in metallicity as a function of doping 

In Kabanov's paper [5], the spectral function a F{o))<xO) in good metals, while in 

poor metals a F{g)) x ft) . In the time evolution of the electronic temperature Te, this effect 



gives rise to the prefactor 



finkl 



in good metals, and 



in 



in poor metals. 



This moderation of the electron-phonon interaction, by the change in metallicity with doping, has 
been accounted for in Equations (1) and (2), by assuming a "two-fluid model", where Cj is the 
linear combination of these two prefactors: 



Ci = (l-a) 



in 



+ a 



n-nkl 



(17) 



with the "good-metal" fraction a given by 



a{p) = a + bp. 



(18) 



We assume that a varies linearly with doping level p. This linearity is justified by the 
approximate linear dependence, with doping, of the residual resistivity of Bi-2212 [6] — see 
Figures 5 and 6. This "good-bad-metal" picture, though phenomenological, is consistent with our 
picture of Bi-2212 becoming more metallic with increasing doping. In the underdoped region, the 
hole carriers are so few that the transport is difficult. At high temperatures, the hole carriers are 
scattered off the random distributed spin moments so strongly that the resistivity is much larger, i.e. 
"bad" metal. In the overdoped region, the effective bandwidth is increased and the spin fluctuations 
are suppressed significantly. As such, the resistivity is decreased, i.e. "good" metal. 

The paremeters a and b are determined by insisting that 



a(0.05) = and a(0.27) = l. 



(19) 



implying that the material is a bad metal at p = 0.05 and a good metal at p = 0.27 . These two 
end-points correspond to the dopings on the superconducting dome when the superconducting 
transition temperature Tc - 0. Figure 1 of Ref. [7] shows that, even at 300 K, the temperature 
where we took our data, these two end-points are still reasonable. 

(2) Doping-dependent penetration depth 

Optical reflectance data of Bi-2212, at 300 K, showed a doping dependence of the 
penetration depth at 12500 cm ', corresponding to the wavelength of the pump and probe pulses 
(800 nm) [6]. 

(3) Doping-independent electronic specific heat 

In Equation (6), the Sommerfeld coefficient y- 8 mJ/mol.K^from Junod [8]. There is no 
doping dependence of y for Bi-2212 at room temperature, as seen from specific heat data of 
Loramef al. [9]. 



50O 



40O 



^ 300 
u 

t 

E 200 

a 



100 



Bj Sr CaCu 




SO 



250 



300 



100 150 200 
Temperature (1^ 

Figure 5. The dc resistivity, p^^(T), extracted 

from low-frequency extrapolations of normal 
state optical conductivities. The dashed curve 
through the data points for the most overdoped 
sample is a fitted curve, p,,,(r) = 0.01557"''' . 

Taken from Hwang [6]. Permission granted from 
lOP Publishing Limited. 
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Figure 6. Resistivity (at 300 K) versus 
doping, obtained from Figure 5. The trend 
is roughly linear, which justifies the use of 
Eqn. (18). 



E. Curve-fitting procedure 

For the results presented in this section the bosonic mode energies are chosen as 
described in our main article: Hot phonon energy [10] = 40 meV, cold phonon energy [11] = 20 
meV while both the hot and cold spin fluctuation energies scale with the centroid position of the 
broad bosonic background [6], with the constraint that the energies are 41 meV at optimal doping. 
We start the fitting by assigning initial values to the parameters (Aesf, -^e-ph, ^st, Tph, /st, /ph)- Based 
on these values the system of five coupled ordinary differential equations is solved and the five 

effective temperatures ( T^ , Tj^^, , T j^ ,T^r ,T^r ) are obtained as functions of time t. We then 

calculate the sum of square errors between the electronic temperature (T^) and the experimental 

data over a range of 7 ps. We then search through the parameter space to find the optimum set of 
parameters (minimum sum of square-errors) using the MATLAB built-in function fminsearch. 

Our fitting results based on six fitting parameters revealed that eventually not all six parameters 
are doping dependent. To investigate this in more detail we begin with an unrestricted fit to 
identify those parameters with unambiguous doping dependence. We found these to be the 
coupling functions lesf and Igph, and that the quality of the fits depends most sensitively on their 
values. 



In the unrestricted fits we further found that the parameters /^f, /ph do not show strong doping 
dependences. If we replace them with constants, and perform restricted fits to the remaining four 
parameters (Ae-st, -^e-ph, 7sf, Tph), we find the fitting errors of each individual sample change the 
least compared to the unrestricted fit across all samples. To determine the optimal values for the 
fractions, we perform a parameter space study. We calculate the fitting error of the 5TM with four 
fitting variables (Xesf, ^e-ph, t,/, Tph) at all possible fraction combinations (f^ufph- 0.01 - 0.96 at 



0.05 intervals). This gives us one error surface for each of the differently doped samples. We then 
average over all of them to obtain the final error surface shown in Figure 7. 

We found two minima in the averaged error surface that lie at almost symmetric locations in the 
parameter space. This is expected since the two bosonic channels are of similar form. The 
channels are distinguishable by their mode energies. At the fraction pairs that produce minima in 
the averaged error surface, all the individual error surfaces have a pronounced minimum as well. 
Thus the choice of fractions at the minimum of the averaged error surface is also a very good 
choice for each sample. 
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Figure 7: Intensity Plot of averaged error surface. The averaged error surface 
features two minima at symmetric locations along the /^f and /ph axis respectively. 
Such symmetry is expected because of the similarity of the bosonic cooling 
channels. A dark-blue coloring denotes fraction combinations that yield very good 
fitting results. The dark-red regions indicate large fitting errors. The locations of 
the minima are (/ph,/sf) = (0.06, 0.46) and (0.46, 0.06). 



From the shape of the error surface two more conclusions can be drawn. First, the fact that the 
two minima are located away from the diagonal suggests that the two bosonic channels play 
essentially different roles. If a single bosonic cooling channel, like in a three temperature model, 
were sufficient, the error surface would feature minima along the diagonal (f^f - /ph) only. We 
found that the channel with a large fraction acts fast and equilibrates with the electronic 
temperature within the first picoseconds (0-3 ps). The channel with a small fraction of hot 
modes controls the cooling process at later times (3-7 ps). Only a combination of both cooling 



channels allows good fitting results of up to 7 ps. Secondly, if the fractions were redundant 
degrees of freedom the error surface would not show sharp localized minima. Instead, a flat 
featureless surface would be expected. Our study thus confirms that both cooling channels should 
consist of one hot bath and one cold bath. 

For the fraction pair (/ph, /sf) = (0.06, 0.46), the doping dependence of the coupling constants is 
shown in Figure 2(d) of the main article. Ig ph exhibits a peak at OPT, while le-sf decreases with 
increasing doping. For the other fraction pair (/ph, /sf) = (0.46, 0.06), not only are the doping 
dependences of Aeph and Aesf reversed, the fits are also not as good as the first case. Since, with 
increasing doping, one is "farther" from the antiferromagnetic region, we should expect a 
decreasing le sf- We thus take the doping dependence of the coupling constants based on the first 
fraction pair, that is, (/ph,/sf) = (0.06, 0.46). 



F. Time evolution of electronic, phononic and spin fluctuation subsystems 

Figure 8 reproduces Figures 2(a)-(c) of the manuscript. After the initial photoexcitation, 
the electrons heat up, rapidly transferring their energy to other subsystems. Due to the small 
number of hot phonons (fph - 0.06), equilibration with electrons (T^^l" = J. ) is reached within the 
first picosecond. The value of le ph determines how fast the hot phonons respond to the electrons, 
and characterizes the sharpness of TifJ^ -peak — from UD to OPT leph increases, resulting in a 
sharper hot-phonon peak; in the OD samples the hot phonons react slower due to a smaller 2e-ph- 
Hot spin- fluctuations (r^-[^ ), on the other hand, equilibrate with electrons within the first 2 ps, due 

to their larger population (/ph = 0.46) compared to hot phonons. The subsequent cooling of the 
electronic (Te) peak is thus mainly determined by the coupling between electrons and hot 
spin-fluctuations. With increasing doping, lesf decreases, which is reflected in a less-sharp T^^^ 

-peak. Between 1.5 to 3 ps, cold spin-fluctuations (r,f;^ ) equilibrate with electrons. Beyond 3~ps 

(hot and cold) spin-fluctuations are not capable to decrease Tg any further — further electronic 
onwards cooling is due to cold phonons only. At times >7 ps all five temperatures eventually 
equilibrate and the tail flattens out. The small shoulder at -0.5 ps is an extrinsic effect present in 
all UD samples, but does not affect the quality of the fits. 



G. Error Analysis 

To estimate the fitting error of our parameters we repeat the same fitting process for a 
large number of experimental data sets. While the actual parameters are obtained by fitting the 
average of the experimental data, the respective error bars are estimated by fitting all individual 
data sets separately. We then took standard deviations of the distributions of the various fitting 
parameters to obtain the error bars. 

Figure 8 shows the original experimental data as grey shaded area behind the averaged curve used 
to obtain the actual fitting parameters. When we fit the wide range of noisy data we observe that 
only a narrow range of fitting parameters is necessary to obtain good fits. This suggests that our 
model parameters react sensitively to changes in the shape of the curve which is an indicator for a 
tightly constrained model. 
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Figure 8: Open circles = averaged data 
multiplied by a proportionality factor. Solid 
lines are time evolutions of the temperatures of 
the various subsystems: electrons (black), hot 
phonons ( T,,",'; , red), cold phonons ( TJ^,';, , 

green), hot spin fluctuations (T^f, , blue), and 

cold spin fluctuations ( T'^,^ , magenta). 

Reproduced from Figures 2(a)-(c) of 
manuscript. Grey shaded area denotes the 
individual experimental data sets. 



H. Other choices of hot phonon and hot spin fluctuation energies 

Besides the scenario presented in the main text, we show here the doping dependences of 
Ae-ph and /le-st for three other combinations of hot phonon and hot spin fluctuation energies (Figure 
9). The fraction pairs used for these scenarios were obtained from error surface studies as well. 

(a) Hot phonon energy = 40 meV [10]. Hot spin fluctuation energies taken from the energies 
of the magnetic resonance mode from neutron scatterig data [12], and the positions of the 
peak in the bosonic spectrum from optical spectroscopy data [6]. 

(b) Hot phonon energy = 70 meV [10]. Hot spin fluctuation energies taken from the energies 
of the magnetic resonance mode from neutron scatterig data [12], and the positions of the 
peak in the bosonic spectrum from optical spectroscopy data [6]. 

(c) Hot phonon energy = 70 meV [10]. Hot spin fluctuation energies taken from the centroid 
of the broad background of the bosonic spectrum in optical spectroscopy data [6]. 



As shown in Figure 9, the doping dependence of le-ph and le-sf is quite robust against 
whichever scenario as listed above is used. 
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Figure 9: Doping dependences of }^e-ph 
and >ue-sf corresponding to scenarios (a) - 
(c). Comparing these three figures with 
Fig. 2(a)-(c) of the main text, we see that 
the doping dependences of le-ph and 2e sf 
are preserved. Only their magnitudes 
have changed. 



The doping dependence of Ae ph and Ig-sf, obtained from our fits looks physically reasonable — 
Ae-ph is roughly BCS-like, peaking at optimally doping, while Ag st gets stronger while approaching 
the antiferromagnetic insulator region. We consider this an indirect proof of our model. If the 
five-temperature model were not vaUd, then there is no reason, a priori, that the doping 
dependence of 2e ph and Ae sf would look the way they appear. In particular, in the underdoped 
region, most of the electrons are localized as local moments, which are responsible dominantly 
for the spin fluctuations, while the doped holes are effectively low-energy charge-carrying 
quasiparticles. The differentiation of the spin and charge degrees of freedom leads to the effective 
quasiparticle-spin fluctuation coupling, in addition to the electron-phonon coupling. More 
importantly, we were encouraged to suggest the five-temperature model by the fact that this new 
model interprets the experimental data nicely throughout the wide doping region considered in 
the present work. 
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